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1 Introduction 


Recently, the study of electroweak processes in ultra-strong magnetic fields became extremely 
actual. There are several cases when such fields could be possible. 

The conception of magnetars [T] (a kind of neutron stars with extremely high magnetic helds) 
encouraged the scientists to reanalyze the electroweak processes in external magnetized media 
with large scale magnetic helds. According to the McGill Online Magnetar Catalog [2], there 
are 28 objects expected to be magnetars. Approximately a half of them are the soft gamma-ray 
repeaters (SGR) and the others are the anomalous X-ray pulsars (AXP). The above-mentioned 
magnetic held scales are greater than the so-called Schwinger held limit—a natural scale for 
the magnetic held strength^ B^, = m^je ~ 4.41 x 10^^ G. Even stronger magnetic helds could 
exist in the early Universe. Note that for the electric held instead of the magnetic held the 
shown limit is crucial because of the intensive vacuum electron-positron pair creation that is 
equal to the short circuit of an electric held generator. On the other hand, the vacuum is stable 
in a magnetic held, so its strength could exceed the critical value Be- Moreover, the magnetic 
held could play a stabilizing role in the case when it is perpendicular to the electric one. In 
such a conhguration, the electric held could be even greater than B^. The invariant form of the 
stability condition is as follows: 

where F^‘' is the electromagnetic held tensor. 

Nowadays, in our Universe from stars to galaxy clusters, magnetic helds of different scales 
from mG to G are observed. There is an open question about the causes of observable large scale 
magnetic helds. On the one hand, the magnetohydrodynamic mechanism explains successfully 
the rate of the observed helds, but for its realizations some seed helds of order 10“^^ -10“^® G 
are required [3]. In this case, the question emerges when and how do the hrst seed magnetic helds 
appear? In the early Universe, in the interval between the stages of the QGD phase transition 
and nucleosynthesis, very strong magnetic helds, the so called “primary” helds, in principle, 
could exist, with an initial strength of the order 10^^ G [1] and even more ( 10^^ G) [5]. Their 

evolution during the expansion of the Universe could determine the existence at the present 
stage of coherent large-scale (~ 100 kpc) magnetic helds with an intensity I'N-' 10-21 G. These 
helds, in turn, could be enhanced by the galactic dynamo mechanism to the observed values 
of galactic magnetic helds ~ 10“® G. Possible origins of primary strong magnetic helds and 
dynamics of their evolution in the expanding Universe are the subject of intense research (see, 
for example, the surveys [6H8] and references cited therein). 

Finally, one more case when strong magnetic helds could be possibly observed are the 
contemporary elementary particle accelerators. The fact of observing of the helds of magnetar 
scales in the conditions of the Earth laboratories seems to be suprising but the helds of such 
magnitudes could really be formed during non-central heavy-ion collisions [S]. 

While solving the series of principal problems of the interactions of charged particles with an 
electromagnetic held, the method has gained a great value, in which the impact of external helds 
is considered on the base of the direct solutions of wave equations in the external electromagnetic 
held instead of using the perturbation theory. The method was hrst proposed by Sokolov and 
later developed by Sokolov and Ternov while solving the problem of the synchrotron radiation. 
The corresponding idea is now known as the Furry picture. The calculation of the specihc 
physical phenomena supposes the using of the Feynman diagram technique with the following 
generalization: in the initial and hnal states, a charged particle is in the external held and 
is described by the solution of the wave equation in this held, and the internal lines of the 
charged particles correspond to the propagators built on a base of that solution. The detailed 
description of the calculation technique of the processes in external helds one can hnd, for 

^ We use the Planck units: ft = 1, c = 1, and the Minkowski metric with the signature (H—-). 
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example, in reviews and book, see e.g. Refs. unHis]. The method is very useful in the case of 
strong magnetic fields when the account of the field influence can not be made with the help 
of the perturbation theory. By virtue of the vacuum stability in superstrong magnetic helds, 
it is possible to study processes in the fields with the intensities significantly above the critical 
value Bg. 

The extent of influence of the external field on the propagation of a charged particle is 
determined by its specific charge i.e. the ratio of the particle charge to its mass. So, the light 
charged fermions especially electrons are mainly sensitive to the external field influence. 

The expression for the exact electron propagator in the constant uniform magnetic field was 
first obtained by J. Schwinger m in the Fock-Schwinger m proper time formalism. There 
exists a number of works where another forms of propagators were derived. For example the 
case of a superstrong magnetic field was analyzed in Ref. [18] and the contribution from the 
ground Landau level to the electron propagator was obtained. In Ref. m, the propagator was 
transformed from the Schwinger form [16] to the expansion over all Landau levels. In Ref. [20] . 
the expansion of the electron propagator in the power series of the intensity of a magnetic field 
was presented. The exact propagator of an electron in a constant uniform magnetic field as the 
sum over Landau levels was obtained in Ref. m from exact solutions of the Dirac equation by 
a direct derivation of the standard methods of the quantum field theory. 

In our opinion, a knowledge of different representations of the charged particle propagators 
in an external magnetic field and the conditions of their applicability would be very important. 
There exist several precedents when misunderstanding of such conditions led to incorrect studies. 
For instance, a calculation of the neutrino self-energy operator in a magnetic field was performed 
in Refs. [221123] by analyzing the one-loop diagram v ^ e~ W~^ —>■ v. The authors restricted 
themselves by the contribution to the electron propagator from the ground Landau level. As 
there was shown in Ref. [23] , in that case the contribution from the ground Landau level did not 
dominate due to the large electron virtuality, and contributions from other levels were of the 
same order. Ignoring such a fact led the authors [221123] to incorrect results. Another example 
of this kind was an attempt to reanalyze the probability of the neutrino decay v ^ e~ W~^ in an 
external magnetic field in the limit of ultra-high neutrino energies, calculated via the imaginary 
part of the one-loop amplitude of the transition v ^ e~ W~^ —>■ v. Initially, the result was 
obtained in Ref. [25], later the other authors [26] repeated the calculation and insisted on 
another result. The third independent calculation m confirmed the result of Ref. [25]. The 
most likely cause of the error in Ref. [26] was that the authors restricted themselves by only 
linear terms in the expansion of the VF-boson propagator over the electromagnetic tensor 
whereas the quadratic terms were essential as well. 

Among the works dedicated to the exploration of the propagators of charged particles in the 
external field one should note Ref. [28], where the calculation of a neutrino self-energy operator 
in a magnetic field in an arbitrary ^-gauge was held. It was demonstrated that in spite of the 
fact that the self-energy operator depended on the gauge parameter the observed neutrino 
properties arising from its dispersion law, as expected, were gauge-invariant. 

Analyses of literature showed that there are no methodically important expressions for 
propagators of charged W- and <l>-bosons as expansions over Landau levels. In the present work 
we calculate, for the sake of completeness, various representations of propagators of charged 
particles in a constant uniform magnetic field, that could be useful for the analyses of the 
electroweak processes in strong magnetic fields. Different methodically important transitions 
from one representation to another are also performed. In Sec. [21 the general forms of the 
electron, charged W- and <h-boson propagators in the Fock-Schwinger proper time formalism 
are presented. In Sec.O several important comments on the translation and gauge noninvariant 
phases are made. Further in Sec. HI expansions are presented of propagators in the weak-field 
approximation. In Sec.jSl the procedure of transformation of the electron propagator to the form 
where it is presented as the sum over all Landau levels is shown. The formulas for propagators 
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of the charged W- and scalar charged <I>-bosons in an arbitrary ^-gauge expanded over Landau 
levels similarly to the electron propagator expansion are derived in Sec. [6] for the first time. 
Some details of calculations are presented in the appendices. 

We use the notation for the 4-vectors = [t,x,y,z)] the index numbers designate the 

different 4-vectors. It is implied in the paper that the squared masses have infinitely small 
imaginary parts m? —>■ rm? — ie. 


2 Propagators in the Fock-Schwinger proper time formalism 

The electron propagator in a constant uniform magnetic field in the Fock-Schwinger proper 
time formalism can be presented in the form 

-X 2 ). (1) 


Here, S{X) is the translation and gauge invariant part of the propagator: 


S{X) = - 


00 

i/3 f ds 


1 

2(47r)2 J ssin(/3s) \ s 
0 


/ 


cos{f3s){X(f(f'y) — isin(/3s)(W(^7)75 




sin(/3s) 

X exp < — i 


{X(p<f-f) -I- rUe 


2 ^{X^^X) 
m„s + 


2cos(/3s) — sin(/3s)(7¥?7) 


4s 4tan(/3s) 


{Xipy^X) 


( 2 ) 


where /3 = eB, e is the elementary charge, me is the electron mass, = {Xi — ^ 2 )^, (pap = 
Fap/B is the dimensionless electromagnetic tensor, pay = ^ is the dual dimensionless 

tensor _ -gi); tensor indexes for 4-vectors and tensors in parentheses are 

assumed to be contracted consecutively, for instance: {Xppj) = X°‘pai3^^^'y^i- 

Integration with respect to the variable s (corresponding to the proper time) in the expres¬ 
sion m needs to be redefined due to the poles of the integrand in the points s = 'nk/fi, where 
/c = 0,1, 2... It is assumed that the integration is performed in a complex plane s over the 
contour that begins in the point s = 0 and lies below the real axis closely to it. This contour 
may also be rotated to the negative imaginary semiaxis (see below). 

The phase 4>(Xi,X2) in the formula ([T]) is translationally and gauge noninvariant and can 
be defined in terms of an integral along an arbitrary contour as: 


X 2 


$(Xi,X2) = -e / dX^X^(X), 

(3) 

J 

Xi 


K^^{X) = A^^{X) + - X2)^ . 

(4) 


Similarly to the equation o, one can define the propagators of the IT-boson and the charged 
scalar <I>-boson in a magnetic field (note that we assume negatively charged W~ and 4>~-bosons 
to be particles): 


G;,^)(Xi,X 2) = e'‘^(^i’^2) G^,{Xi - X 2 ), (5) 

D('^)(Xi,X 2 ) = e^^(^i’^2)7)(Xi -X 2 ), (6) 

where the phase 4>(Xi,X2) is defined by the same formulas (j3|l, (|4|l. 


4 













It is often convenient to use the Fourier transforms of the translationally invariant parts of 
propagators: 


S{X) = 
G^,{X) = 
D{X) = 


(27r)4 


5(g) 


(7) 

( 8 ) 
(9) 


One can obtain the Fourier transform of the electron propagator from Eqs. ([2|) and d?]) in the 
form: 


OC 

S{q) = [■ 

J ' 


ds 


cos{f3s] 


exp 


2 2,2 tan(/3s) 


[(97)11 + 


rrip. 


cosips) - sin(/3s) 


( 97)1 
cos(/3s) J ’ 


( 10 ) 


where the following denotations for multiplications of longitudinal and transverse 4-vectors are 
used (g 7 )|| = (qipif'y), (g7)± = (g 99757 ). In the coordinate system where the 3d space axis is 
directed along the vector of an external magnetic field B, the 4-vectors with the indexes _L and 
II belong to the Euclidean {1, 2}-subspace and the Minkowski {0,3}-subspace respectively. Eor 
instance, p± = ( 0 ,pi,p 2 , 0 ), p\\ = (po, 0 , 0 ,p 3 ). 

The Eourier transforms of the propagators of the VE-boson ([5]), ([8]) and the charged scalar 
<I)-boson ([6]), ([9]) depend on the gauge selection. In an arbitrary ^-gauge they take the form 


G^^uiq) = - 


tan(/3s)/(/3s)) 


cos(/3s) 


(1 - cos(2/3s)) 


sin(2,ds) 


m 


W L 


9 /i + (‘Pq)/! tan(/3s) ] { qu + {qp>)u tan(/3s) 


+i f tan(/3s)) (e'^^w _ | , 

00 


_ds_^-is(5m^-g|+g2 tan(/3s)/(/3s)) 

cos(/3s) 


( 11 ) 


( 12 ) 


In the Eeynman gauge, where ^ = 1, the Eourier transform of the kE-boson propagator signifi¬ 
cantly simplifies |29j : 

00 

J cos(/3s) 

0 


Qfiu + {p>p>)fiu (1 - cos(2/3s)) - sin(2/3s) 


(13) 


The Eourier transform of the propagator of the charged scalar <I>-boson in the Eeynman gauge 
has the form: 


Diq) = 


ds 


cos(/3s| 


exp 


2 2 , tan(/3s) 2 

—IS ( rriy^/ — gy H —— g_|_ 


/3s 


(14) 
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3 Comment on the noninvariant phase 

The expressions ([3]), (H]) for translationally and gauge noninvariant phase ^{Xi,X 2 ) written in 
a covariant form look rather cumbersome. Some authors prefer to write down the phase in more 
compact form, fixing the gauge by choosing 4-potential as A^^{X) = {0,0,x B,0), to obtain: 

pB 

^X,X') = - — {x + x'){y-y'). (15) 

Nevertheless, the exactly covariant form of the phase ([3]), dH) appears to be much more 
efficient in an analysis of the closed loops containing several propagators of charged particles. 
Note that by virtue of the property of the 4-vector ([3]): d^K^ — duK^ = 0, the path of integration 
from Xi to X 2 in Eq. ([3]) is arbitrary. In the case of the two-vertex loop, the sum of the phases 
contributed to the amplitude is zero: 


$(Xl,X2) + $(X2,Xi) =0. (16) 

In the presence of the three or more vertices, there appears a nonzero total phase from all 
propagators, which is nonetheless translationally and gauge invariant. This can be easily shown 
by rewriting the 4-potential of the constant uniform external field in an arbitrary gauge in the 
following form: 


A^(X) = ^X,F"^ + 9^X(^), (17) 

where x(^) is an arbitrary function. Given Eq. (1171) . one automatically obtains that d^A'^ — 
d'^A^ = . Integrating ([3]) with taking into account the relation (IIZD, we obtain 

HXi,X2) = - | (X1EX2) - e [x{X2) - X(^i)] . ( 18 ) 

It is seen from Eq. (jlSp that during the summation of phases in a closed loop, the terms 
containing the function x is completely cancelled, providing the gauge invariance. It can be 
easily checked that the sum of the expressions (fTSl) in a closed loop would be translationally 
invariant also. 

For instance, for the three and four propagators in a loop, the sum of phases takes the form: 
$(Xi, X2) + $(X2, X3) + 4>(X3, Xi) = (Xi - X2)^F^"(X2 - X3). , (19) 


$(Xi,X2) + ^(X2,X3) + $(X3,X4) + 4>(X4 ,Xi) 

= -|(Xi - X3)^F'^"(X2 - X 4 ). . ( 20 ) 

In a general case, the sum of n phases can be presented in the following form: 


^tot = -lY.^X,FX,+,) 

i=l 


X„+i=Xi 


e 

2 


n—11—1 




where 


Zi — Xi — Xj+i. 


( 21 ) 
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4 Propagators expanded over the weak field 


Using the exact expressions m and (I13p makes the calculations sufficiently cumbersome. At 
the same time, magnetic fields existing in Nature, except the early Universe, are always weak 
compared to the critical field for the kU-boson, Bw = triple ~ 10^“^ G. Hence propagators of 
the kU-boson and the charged scalar <l>-boson could be expanded in the series over the powers 
of /? = eB assuming /3 to be a small parameter. Keeping the terms up to the second order for 
the kk-boson propagator the in Feynman gauge one can get: 


■ ,2 ^ 2 ^, (,2 _ „ 2^)2 

2 

+ — -Vtt 1 + 4 


2 


+ i/32 




(g2 - mlrY 


1 


(g 2 - m 2 ) 3 _ 


+ Oif3^ 


( 22 ) 


It is easy to perform the similar calculation for the kU-boson in an arbitrary ^-gauge, but the 
hnal expression is too cumbersome. 

Comparing (fT^ and (fT^ one can see that the <I>-boson propagator D{q) differs only by the 
sign from the coefficient before in the expansion of the propagator G^y{q) over the three 
independent tensor structures. So, 


D{q) = 


q^ - mfy 


-i/3^ 


(q2 — m,2 


+ 


2 ql 




(g2 — m,2 




+ 0(13^ 


(23) 


The asymptotic expression for the electron propagator S{q) in the case when the field 
intensity is the minor dimensional parameter of the problem /3 <C m^ <C m^ can be obtained 
in the similar way. In the case of the weak field approximation, the expansion for the electron 
propagator can be derived in the form m-- 


Siq) 

+(3^ 


^ . (gy) + me ^ ^ (g7)|| + "i. 


2 i 


2(,2_„2)2 

(^n - 'ml) (g 7 )± - qI (( 97)11 + 'me) 


g2 — ml 


(g2 - m2)4 


+ 0(/3^). 


(24) 


Note that in this expansion, the contribution from the region of the small virtual momenta, 
g2 ~ ml <C increases in each next terms. Wherein in the case of the moderately strong 
field, ml P m^y, the expansion ([24]) is inapplicable and it is necessary to use the exact 
expression (fTn|) for the propagator. 


5 Electron propagator as an expansion over the Landau levels 

Let us rewrite the Fourier transform m of the translationally and gauge invariant part of the 
propagator introducing a new variable of integration v = f3s: 

00 

Siq) = ^ J dv exp(-ipu) | [(g 7 )|| + me] fi{v) 

0 

- [(97)11 + me] f 2 {v) - (g 7 )x /3(u)| , (25) 

where we use the notations: 


7 


















( 26 ) 


fi{v) = exp(—ia tanu), 
f 2 iv) = tant; exp(—ia tanu), 

fsiv) = - 0 — exp(—ia tant>), 

and also p = {ml — gy)//3, a = q\/fi- Because of the periodic character of the functions 
(j = 1 ) 2 ,3), fj{v) = fj{v + nvr), we can slice the integration area ( 0 , oo) to the intercepts 
(0, tt), (vr, 27r),... (utt, (n + l) 7 r).... Making in each intercept the variable substitution v —>■ 
V + nvr, we can rewrite: 


n =0 


du exp(—ipu) fj{v) = exp(—ipn-Tr) j du exp(—ipu) fj{v) 
1 

1 — exp(—ipvr) 


A 


■3 ’ 


(27) 


where 




TT 

J du exp(-ipu)/j (u). 


(28) 


It suffices to compute the integral Ai because the other two integrals can then be found 
with the help of the expressions: 


. ■ 9 A 

M = 1 — , 

oa 

Tls = -- (l-e-‘^")-^^i. 
a a 

One can easily verify the validity of the last formula substituting A 3 in the form: 


(29) 


A?, = 


TT 

— J dv exp(—ipu) — ^exp(—iatanu)^ 


(30) 


and computing the integral by parts. 

It is not difficult to calculate the integral Ai (see Appendix A for details). Appling the 
formulas ([25l), ([26]), ([28]), ([29]) and the exact expression for Ai (06]), we finally write the Fourier 
transform of the traslationally and gauge invariant part of the electron propagator in the form: 


sii) = E 


^ - m2 - 2n/3 


[(97)11 + 


m. 


dn{a) - X ( 77 ^ 7 ) (a) 


- (97)± 2n 


dn{a) 


a 


where a = q\ll5 and 


dn{a) = (-l)”e-“[L„( 2 «) -L„_i( 2 a)] 


(31) 

(32) 
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6 Propagators of the charged W and ^-bosons expanded over 
the Landau levels 

In a case when a magnetic field is sufficiently strong, B > Be = m\le^ it is convenient to use 
the electron propagator expanded over the Landau levels. The procedure for obtaining this 
representation is described in Ref. |19] (see also Refs. [201121] ). Note that there is a misprint in 
the expression for the propagator obtained in Ref. [19] . namely the term on the second line in 
formula (4.33) must contain the extra factor (—i). This misprint was corrected in later works, 
Ref. [20] (formulas (39) and (40)) and also Ref. [30| (formulas (13) and (14)) but without any 
comments. 

The propagators of the W and <h-bosons can be similarly expanded over the Landau levels. 
As it was noted in the Introduction, in the conditions of the early Universe magnetic fields 
could possibly exist with the scales of the critical magnetic field value for the lU-boson, Byi/ = 
rn^/e ~ 10^“^ G. In this case, a knowledge of the vector boson propagator in the form of 
expansion over the Landau levels may be useful. 

The Fourier transform of the traslationally invariant part of the VF-boson propagator in an 
arbitrary ^-gauge is shown in Eq. (HH). Let us rewrite Eq. m in a more convenient way for 
further computations: 


Gi,u{q) = -^ j dve ^ 


—ipv 


fi{v) - f 2 {v) - fsiv) 


13 J 


+ ) fi{v) 


(33) 


+ { + q^,{q(p)u - i ^ ] fiiv) + {q>q)fM{q^)v h{v) 


where we introduce the functions: 


h{v) 

Mv) 

h{v) 

h{v) 

h{v) 


1 


cos V 
cos(2 v) 
cos V 
sin(2 v) 
cos V 
tanu 


exp(—ia tanu), 
exp(—ia tanu), 
exp(—ia tanu), 


d 


exp(-ia tanu) = i— fi{v ), 
cos V oa 

tan^ V , . , / N 

-exp(—la tanu) = — Jiiv) 

cos V oa^ 


(34) 


and also designate p = (m^ - g|)//3, pg = (C rn^ - g|)//3, a = q\/l3. 

Taking into account that fj{v + vrn) = (—l)"'/j(u) (j = 1,2,3,4,5), and slicing the inte¬ 
gration area (0, oo) to the intercepts (0, vr), (vr, 27r),... (nvr, (n -|- l)7 r)... and making in each 
intercept the variable substitution v ^ v + nvr, we can rewrite: 


OO 

j du exp{-ipv) fj{v) 
0 


1 

1 + exp(—ipvr) 


A 


J ’ 


(35) 
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where the integrals were introduced: 




TT 

J dv exp{-ipv) fj{v), (j = 1,2, 3,4, 5). 


( 36 ) 


In order to obtain the final expressions for Ai, A 2 , A^, it is necessary to use the relation 
Yl^=o^n{cAj = 1 (assuming that L_i(x) = 0). To present the results in a more compact way 
one should employ the notations: 


2 a 


The integrals take the following forms (see Appendix B for details): 

(o) 

p + 2n — 1 ’ 


ii =-2i(l + e-^^^) ^ 


(37) 


n=0 


A 2 = -i (1 + ^ ^n- 2 (a) 


n=0 


p + 2n — 1 ’ 


A.s = — 




n=0 


^n(^) ^ 71 — 2 ( 0 ^) 

p + 2n — 1 


00 Q! 


A4 = 2 (1 + e-'^^) 


n=0 


4-i(«) 

p + 2n — 1 ’ 


00 «// 


is =2i(l + e-^")^ 


n=0 


C-i(«) 

p + 2n — 1 


Substituting the values of ii, ...is into the expression for a propagator ([331) we finally 
obtain: 


G^ig) = E 


— 1 


^^gj-m^-/3(2n-l) 


2 (p 'p )1 (o:) 


(ptp)^^ ^7^,(0:) + £n— 2 {oi)^ T i :p^j/(o:) ^n—2(0:)^ 

2 QfiQi/ T i/3 ^n—iicn) 


qfl -irn^ - /3(2n - 1) 

+ i( 2(pg)^gfi. + 2q^{qp)y - i/3 (pp);,i. ) ^n_i(a) 


-2{pq)^{qp)i,C_4{a) 


(38) 


Note that unlike to the electron propagator, the contribution from the ground level n = 0 
into the propagator of the kk-boson has a special feature. This contribution has the gauge 


10 












independent form: 


OfM) - A ^ a 0-'''^'“ ■ ( 39 ) 

y|l ^ 

It is essential that it contains a pole in the point — j5. So, if the magnetic field rises 

up to the critical value for the Vh-boson, Bw = rn^/e ~ 10^^ G, then the so-called instability 
arises of the perturbation theory for a W boson vacuum, see Refs. [3THM1- 

The propagator of the <l>-boson in an arbitrary ^-gauge, D{q), can be obtained from Eq. ([38|1 . 
similarly to the case of the weak field approximation, in the form: 


D{q) = 


2iin-i{a) 


n=0 


(2 ^ - 1) 


(40) 


Note that the summation over n in the formula (140 p starts formally from n = 0, but actually 
from n = 1, because ^-l{a) = 0 by definition. This means that the propagator of the boson, 
as one could expect, does not contain a pole at the point ~ 


7 Conclusion 

Different representations of the propagators of charged particles in a constant uniform magnetic 
field have been analyzed. The expressions for the propagators of the vector kT-boson and the 
scalar <I>-boson of the Standard Model in an arbitrary gauge expanded over the Landau levels, 
Eqs. (I38P and (|i0]l . were derived for the first time. 
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Note Added in Proof 

After the paper was completed, we have been aware of the paper |35] where the formula for 
the propagator of a charged scalar particle expanded over Landau levels was obtained. Our 
expression for this propagator, taken in the Feynman gauge, coincides with the result of Ref. [35] . 
to the notations. We thank A. Ayala for informing us on their paper. 


A Calculation of the integral Ai 

To compute Ai let us rewrite fi{v) in the following form: 


fiiv) 


exp(—ia tanu) = exp 



-e-2i^ + 1\ 

_e-2ii; _ 1 J 


(41) 


The right-hand side of this equation can be expressed through the Laguerre polynomials: 


Ln{x) 


1 X 

— 6"^ - 

n! dx”' 


(x' 



(42) 
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The generating function for the Laguerre polynomials is determined by 


1-t 


exp 


X t 

1-t 




n=0 


for hi < 1, so 


exp 


/ f \ °° 

(-T^ 

^ ^ n=0 


( 43 ) 


(44) 


where we assume L_i(x) = 0. Denoting in the right part of the formula (j4ip : —e = t, 
2a = X, and using the identity 


exp o 


X t+ 1 


2 t-l 


= exp 


xt \ f X 
exp -- 


1 -t 


2 / ’ 


(45) 


let us transform Ai to the form 


^ oo 

Ai= due“" ^ 

^_ r\ 


71=0 *- 


Ln{2a) - Ln-i{2a) 


(—I)”" exp(—2inu) exp(—ipu) 


= e 


^ (-!)"■ [L„( 2 q:) - L„_i(2a)] J du exjp[-i{p + 2n)v] 

0 

(-1)' 


n=0 


= h _ 




71=0 


p + 2n 


Ln(2a) - Ln-i{2a) 


(46) 


B Calculation of the functions Aj 

In order to calculate the functions Aj, it is worthwhile to introduce the auxiliary integrals: 


C(a) 

= J exp(—i/9u) 

0 

exp( 

S(a) 

TT 

= J exp(—i/9u) 

0 

exp( 

E(±)(a) 

= C(a)±iS(a) = 

TT 

1 du 


cos V, 


smu, 


Ai can be represented in the form 

TT 

— / du exp(—ipu) cosu — (exp(—iatanu)) 
a J av \ / 


Ai = 


(47) 


and further, integrating by parts we obtain: 


Ai = 


a 


— 1 — exp(—i/?7r) + \pC{a) + S(a) 


(48) 
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The integrals A 2 and can be expressed through Ai, C(q;) and S(a): 


A 2 — 2 C(q;) — Ai, 


( 49 ) 


i3 = 2S(a). (50) 

In order to calculate C(a) and S(q;), let us compute £^^^( 0 ;) and use the expressions: 


C(a) = i E^+\a) + E^-\a) 


S(a) = ^ E^+\a) - E^-\a) 


The integral (a) equals to 


E^^^(q;) = —i [1 + exp(—ipvr)] 


dn (q:) 


n=0 


The functions dn{a) are determined by the equation 
The integrals C(q;) and S(a) can be expressed as 


. LXJ 

C(a) = - ^ [1 + exp(-ip7r)] ^ 


n=0 


p + 2n =F 1 


dnipd) T dYi—\{ot) 
p + 2n — 1 


(51) 

(52) 

(53) 


(54) 


Qt \ 1 n I ( ■ dn{a) — dn-i{a) 

S(a) = - ^ [1 + exp(-ip7r)] ^ - 


n=0 


p + 2n — 1 


(55) 
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